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Abstract. The most frequently used in physical application diffusive (based on the Fokker-Planck equation) 
model leans upon the assumption of small jumps of a macroscopic variable for each given realization of 
the stochastic process. This imposes restrictions on the description of the phase transition problem where 
the system is to overcome some finite potential barrier, or systems with finite size where the fluctuations 
are comparable with the size of a system. We suggest a complementary stochastic description of physical 
systems based on the mathematical stochastic storage model with basic notions of random input and 
output into a system. It reproduces statistical distributions typical for noise-induced phase transitions (e.g. 
Verhulst model) for the simplest (up to linear) forms of the escape function. We consider a generalization of 
the stochastic model based on the series development of the kinetic potential. On the contrast to Gaussian 
processes in which the development in series over a small parameter characterizing the jump value is 
assumed [Stratonovich R.L., Nonlinear Nonequilibrium Thermodynamics, Springer Series in Synergetics, 
vol.59, Springer Verlag, 1994], we propose a series expansion directly suitable for storage models and 
introduce the kinetic potential generalizing them. 

PACS. 05.70.Ln Nonequilibrium and irreversible thermodynamics - 05. 40. -a Fluctuation phenomena, 
random processes, noise, and Brownian motion - 05.10.Gg Stochastic analysis methods (Fokker-Planck, 
Langevin etc) 



1 Introduction 

One of the aspects of modelling the behaviour of a complex 
physical system consists in introducing a random process 
capable of describing its essential properties. The most 
common (and in practice almost unique) class of stochas- 

' tic processes where reliable results can be obtained is the 
class of Markov processes. Said processes in their turn can 
be subdivided into different families. The most widespread 

' is the model of diffusion process with Gaussian noise su- 
perimposing on the macroscopic dynamics. The Poisson 
random processes (or " shot noise" ) present a second bench 
point [T] together with the former subclass covering the 
most common physical situations. In the present paper 
we bring into consideration the stochastic storage models 
based on essentially non- Gaussian noise and treat them 
as a complementary alternative to the diffusion approxi- 
mation (that is to the Gaussian white noise) . We consider 
the phase transitions in such models which resemble the 
noise- induced phase transitions pQ. 

The class of stochastic storage models 2,3^ presents 
a rather developed area of the stochastic theory. As op- 
posed to the common diffusion model [HE] it contains as 
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essential part such physical prerequisites: (i) limitation of 
the positive semispace of states, (ii) the jumps of a ran- 
dom physical process which need not be considered small; 
(iii) essentially non-zero thermodynamic flux explicitely 
specified by the process of random input. 

The following material allows to conjecture that said 
models provide more handy facilities of describing the 
noise-induced phase transitions than diffusion ones [1]. 
One of the reasons in favour of this is the fact that typical 
probability distributions there are not Gaussian but rather 
exponential and gamma-distributions which are charac- 
teristic for, e.g., Tsallis statistics. The approach of the 
present work is invoked to extend the range of applicabil- 
ity of such models. Attempts have already been made to 
apply them to the kinetics of aerosol coagulation [5], to 
the problems of probabilistic safety assessment method- 
ology [5] etc, and to relate these processes to the Gibbs 
statistics and general theory of dynamic systems .7]. One 
more possible application of the storage models consists 
in the possibility of naturally introducing the concept of 
the lifetime of a system (the random time of existence of 
a given hierarchical level) 2,7,8,9 . It was shown 7 j that 
the ambiguity of macroscopic behaviour of a complex sys- 
tem and the existence of concurring evolution branches 
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can be in principle related to the finiteness or infiniteness 
of its average lifetime. 

It is worth mentioning now that (at least) the simplest 
cases of storage models do not require special probabilis- 
tic techniques, and corresponding kinetic equations are 
treatable by means of the Laplace (or Fourier) transform. 
Up to now such models have not gained much recogni- 
tion in physical problems. We believe them to be rather 
promising especially in the approaches based on modelling 
the kinetics of an open system, where the input and re- 
lease rates could be set from the physical background. In 
the present work we do not intend to cover the variety 
of physical situations akin to the storage models. Having 
discussed the form of the stationary distributions for a set 
of input and release functions (Sect. 2) and their relation 
to noise-induced phase transitions we reconsider the for- 
malism of the kinetic potential and fluctuation-dissipation 
relations (FDR) (Sect. 3) and then pass to the problem of 
reconstructing the underlying stochastic process from the 
available macroscopic data (Sect. 4). The material of Sect. 4 
also considers the possibility of generalizing the classical 
storage schemes to cover more realistic physical situations. 
The concluding Section gives an example of an application 
in the context of a practical problem of modelling a nu- 
clear fission process. 



2 Storage model as prototype to phase 
transition class models 

Stochastic storage models (dam models) belong to a class 
of models well known in the mathematical theory of ran- 
dom processes 0E]. They bear a close relation to the 
queuing theory and continous-time random walk (CTRW) 
schemes [lO] . The visualization object for understanding 
the physical ground of such a model is a reservoir (wa- 
ter pool), the water supply to which is performed in a 
random fashion. The random value X(t) describing the 
bulk amount in a storage is controlled by the stochastic 
equation: 



X(t) = X(0) + A(t) - / r x [X{u)\ du 



(1) 



Here A(t) is the (random) input function; r(X) is the func- 
tion of output (release rate). Usually deterministic func- 
tions r are considered. In the simplest case it is constant: 



r x (X) 



(2) 



' a,X(t) > 
0,X(t) = 0. 

The storage model ([I]) is defined over non-negative values 
X > 0, and the output from an empty system is set to be 
zero ((2]). Therefore the release rate from (JT]) is written as 
a discontinuous function 



r x(i) = r (i) - r ( Q +)Xq, 



Xq 



q = 0. 
q>0 



(3) 



(4) 



More complicated input functions can be brought into 
consideration. Analytical solutions are easy to find for es- 
cape rates up to linear [2l[3] 



r{X)=bX; r(X)=a + bX. (5) 

As to the random process A(t) describing the input 
into the system, it can be specified within various classes of 
processes. For our purposes a partial case will be of special 
interest, namely, that of Levy processes with independent 
increments [2lll0j. It can be completely described by its 
Laplace transform: 

E(exp(-eA(t)) = exp(-tcp(e)) , (6) 

where E(. . .) means averaging. The function <p(9) is ex- 
pressed as 



<p(0) = J (1 - exp(-to)) Xb(x)dx (7) 
o 



with 



A = (f(oo) < oo ; p — A / xb(x)dx = <p'(0) ; 



xb(x)dx 



ip(oo) 



(8) 



The function b{x) and parameters (JSJ) have a transpar- 
ent physical meaning clear from the visualized water pool 
picture of the model. Namely, A describes the intensity of 
Poisson random jumps (time moments when there is some 
input into the pool), and b{x) is the distribution function 
(scaled to unity) of the water amount per one jump with 
average value fi . Thus, vdx = Xb(x)dx is the probabil- 
ity distribution of a generalized Poisson process [10] (for 
a "pure" Poisson process there were b(x) = S(x — p^ 1 ))- 
For illustrative purposes the typical choice 



b(x) = fj,e-^ x 



(9) 



will be considered. In this case the function ([7]) has the 
form 



<p(6) 



xe 

n + e 



(10) 



The parameter p gives the average rate of input into a 
system representing thus an essentially non-zero thermo- 
dynamic flow. The basic property of the stochastic process 
under consideration is thus the violation of the detailed 
balance (absence of the symmetry of the left- and right- 
wards jumps). This intrinsic characteristics makes them 
a candidate for the systems essentially deviating from the 
equilibrium (locating beyond the " thermodynamic branch" 
PQ). From this point of view the thermodynamic equilib- 
rium of a storage model is achieved only in the degenerate 
case A = 0, that is for a system which occupies only the 
state X — (of course, the equilibrium heat fluctuations 
are thus neglected). 

Another property of the model consists in the finite- 
ness of jumps, on the contrast to the custom scheme of 
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Gaussian Markov processes with continuous trajectories. 
Therefore such models can be believed to be more ade- 
quate in describing the systems with fluctuations which 
can no more be considered small (for example, systems 
of small size). We recover however the continuous- walk 
scheme setting A oo, p — > oo and keeping p = X/fj, 
finite. In this case the input is performed with an infinite 
intensity of jumps of infinitely small size, that is the sys- 
tem is driven by a Wiener-like noise process with positive 
increments; if we limit the release rate r{X) with linear 
terms (|S|) the process for the random variable then turns 
to that of Ornstein-Uhlenbeck p ]llipi2] , and the storage 
model presents its natural generalization. More specifi- 
cally, one can introduce the smallness parameter /3 _1 [in 
equilibrium situations with Gaussian noise it would equal 
to kT; generally (3 accounts for the environment and noise 
levels in a system and can be related to the parameter 
in the stationary distribution oj st (X) ~ exp(—pUp(X))] 
such that 

\p = \0, b p = 0b(Px), <p fS {B)=Ptp(6/0) i (11) 

from where the Gaussian case is recovered assuming f3 — > 
oo; the exponent in the characteristic function ([6]) acquires 
now the form ip(6) = p8 + 9 2 a 2 / (2/3) of the Gaussian pro- 
cesses with drift. 

It is instructive from the very outset to trace the re- 
lation of the present models to the stochastic noise intro- 
duced by the Levy flights as well as to the processes en- 
countered in the CTRW. The non-Gaussian stable laws are 
described by means of the characteristic function of their 
transition probabilities in the form [9, 10, 13, 14 . T5 | fT6| J T7T . 

exp(— tD\k\ a ) with the Levy index a, the case a = 2 recov- 
ering the Gaussian law. The generalized central limit theo- 
rem states that the sum of independent random variables 
with equal distributions converges to a stable law with 
some value of a depending on the asymptotic behaviour of 
the individual probability distributions 10, 18j. In the case 
of the storage model the characteristic function tp(0) from 
([7]) [where instead of k there enters 9 after an appropri- 
ate analytical continuation in the complex plane] comes in 
place of D\k\ a . The finiteness of y(oc) = A indicates that 
the trajectories of the storage process are discontinuous in 
time. It is understandable that if the functions b(x) from 
([JJH]) have a finite dispersion, the sum of many storage 
jumps will converge to the Gaussian law with a — 2. From 
the physical picture of the dam model, as well as from the 
analytical expressions like ((6]) we can see that the storage 
models present a class of models where a mimic of the 
long-range flights is effectively introduced, likewise in the 
Levy flights, but the nonlocality is achieved by virtue of 
the finiteness of allowed jumps. Indeed, the trajectories of 
the centered process A(t) — pi present a saw-like lines with 
irregular distribution of the jump sizes. Only in the limit of 
big times and scales it can be viewed as a Wiener process 
with the variance (x 2 ) = a 2 t with a 2 = <p"(0). On shorter 
time scales the behaviour of the process models the fea- 
tures akin to the superdiffusion (to the positive semiaxis), 



and, on the contrary, a completely degenerated " subdiffu- 
sion" to the left since the jumps to the left are forbidden. 
In this context the function ip(0) presents an effectively 
varying Levy index which ranges from the superdiffusive 
region (0 < a < 2) to negative meaning suppression of the 
diffusion. The variable Levy index for diffusion processes is 
encountered in the models of distributed-order fractional 
diffusion equations (see. e.g. [TB]). Actually, it is possible 
to bring into consideration from the outset the input fac- 
tions b(x) pertaining to the "basins of attraction" of other 
stable distributions with Levy indices a =fi 2. The storage 
schemes in which the functions b(x) themselves are stable 
Levy distributions with power- like assymptotics |a;|^ Q_1 
are considered in [3]. 

The CTRW processes [T0l[T3"l[17] are characterized by 
the joint distribution of the waiting times and jumps of 
the variable. The stochastic noise in the storage models is 
a narrow subclass of CTRW where the waiting times and 
jump distributions are factorized, and the waiting time 
distribution is taken in a single possible form ensuring the 
Markov character of the process [18] . This suggests a sim- 
ple generalization to a non-Markovian case. Namely, as- 
suming the storage input moments to be distributed by 
an arbitrary law q(t) instead of used q(t) ~ exp(— Xt) 
we arrive at generalized CTRW schemes yielding semi- 
Markovian processes which can be applied for introducing 
the memory effects into a system. 

The solution to the models (jT][6]) can be found ei- 
ther with the sophisticated apparatus of the mathematical 
storage theory [2][3] or directly by solving the appropriate 
kinetic equation (see Sect. 3, 5). A considerable simplifica- 
tion in the latter case is achieved in the Fourier space 
where up-to linear release rates yield differential equations 
of the first order (the situation is similar to the systems 
with Levy flights which are usually treated in the Fourier 
space; note also the analogy to the method of the "Poisson 
representation" in the chemical reactions problems [12]). 

For the constant escape rate ([2]) all characteristics of 
the time evolution of the model are obtained in the closed 
form [2]. We mention just for reference for r x = 1: 



J exp(-st)E {exp(-6X(t))\X ) alt = 
o 

[exp(-flX ) - 9eM-Xov(s))/v(s)} 
[a-9 + <p(9)\ 

where Xq = X(t = 0), <p(0) the same as in ([7]) and r](s) 
satisfies a functional equation 

r){s) = s + ip [<q(s)] , 77(00) = 00. (12) 

However the feature of interest now is the stationary 
behaviour of the models of the class |T]). Even for con- 
tinuous functions r and b(x) the stationary distributions 
ui s t(X) besides the continuous part g(X) can have an atom 
at zero, that is 

u st {X) = P 6(X) + (1 - P )g(X) , (13) 
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where g(X) is a probability distribution scaled to 1, and 
Pq = limt^,ooP{X(t) — 0). The integral equation for g(X) 
from [3] reads as: 



r(X)g(X)=P v(X,oc) 



x 



v{X - y, oo)g(y)dy (14) 



with the measure v(x, oo) = A f b(y)dy. For the expo- 
nential shape of the input function |9]) for which v(x, oo) = 
Aexp(— fix) the equation ()14|1 can be solved for arbitrary 
release functions r(X) > (C is found from the normal- 
ization condition): 



9(X) 



Cexp(-ixX + \J^ 

AX) ' 



(15) 



The condition of the existence of the stationary distribu- 
tions for arbitrary input and release functions from [3] is 
the existence of some wq such that 



oo w+y 

„ up 7 r du 

w>w J J r(u) 

y—Q u—w 



v(dy) < 1 



(16) 



Similarly the expression for Pq can be written for the gen- 
eral case [3]. There is a simple relation 



Po = 



1+ / (r(x))v(dx) 
o 



between the weight of the zero atom Pq and the average 
lifetime (r(x)) (averaged random time of attaining the 
zero level starting from a point x). The presence of the 
non-zero Pq indicates at the existence of idle periods where 
no elements are present in a system. Such periods can be 
characteristic for systems of the small size (in which the 
values of fluctuations of a macrovariable are comparable 
to their averages) [8] and must influence essentially the 
statistical properties of a system, for example, they impose 
limitations on the maximal correlation time. 

The behaviour of the models |T]) admits a pronounced 
property of noncquilibrium phase transitions (change in 
the character of the stationary distribution) which occur 
when one increases the value of the average thermody- 
namic flow (parameter A). The phase transition points can 
be explored by investigating the extrema of the station- 
ary distribution (cf. the analysis of noise-induced phase 
transitions in p]), that is for the case of (fl"5j) - from the 
condition fir(X) — A — dr(X)/dX. For example, for the 
model with constant escape rate r = a we get two types of 
solutions: converging solution for small input rates and the 
pool overflow (no stationary solution exists) if the average 
input per time unit exceeds the output rate. The criterium 
for the phase transition (I16| in this case reads simply as 
p = a. If p > a, no stationary distributions are possible. 
For p < a the stationary distribution possesses addition- 
ally an atom 8 (X) at X = with the weight Pq = 1 — p/a. 
Explicitly for © and r = 1 the stationary distribution 
u>st(X) for p = X/fi < 1 is: 



co st (X) = P a S(X) + (l-PQ){p-X)e-^-^ x , Pq = 1-p. 

(17) 

Consider now the exit function r(X) — bX, b > and 
the input in the form ((9]). This storage system does not 
have an atom at zero, and the stationary probability dis- 
tribution exists for all input rates - there is no overflow 
in the system: 

io st {X) = fi^X^- 1 exp(-fiX)/r(\/b) (18) 

(r(X/b) is gamma- function) . The phase transition is the 
modal change of the distribution function, which occurs at 
X = b, where the distribution changes its character from 
the exponential (~ exp(— X)) to Gauss-like with a maxi- 
mum at X > (Figure 1). This peculiarity of the station- 
ary distribution can be interpreted as a non-equilibrium 
phase transition induced by external fluctuations. Such 
transitions are typical [Hl21| for the multiplicative type 
of noise. They do not have their deterministic analogue 
and are entirely conditioned by the external noise. The 
phase transition at A = b manifests in the emerging of the 
nonzero maximum of the distribution function although 
all momenta of the distribution change continuously. As 
in the Verhulst model PQ the phase transition at A = b 
coincides with a point in which [D(X)} 1 / 2 — E(X), where 
D(X) is the dispersion, E(X) is the first moment of the 
distribution. With the choice of the input function in the 
form b{x) = 4fi 2 x exp(— 2px) the stationary distribution 



co st (X) = exp(-A/6)(-A/&)( 1 - A / fc )/ 2 x 
X (a/6-i)/2 Jv6 _ i[2v /Zax7& ] 

(J is Bessel function). The behaviour of the distribution is 
qualitatively the same as on Fig.l with the phase transi- 
tion point at A = b as well. In both cases we have the phase 
transitions which are caused by an additive noise (which 
does not depend on the system variable). The existence of 
phase transitions for the additive noise is closely related 
to the long-range character of the distribution function 
of the noise and such transitions were discovered in sys- 
tems with, e.g. Levy type of additive noise [91114] where 
the structural noise-induced phase transitions are condi- 
tioned by the trade-off between the long-range character 
of the flights and the relaxation processes in the model. 
Analogous conclusions for other types of superdiffusive 
noises are also drawn in [2H1I1ZT] etc. We can thus state that 
the effective long-rangeness in the storage models leads to 
similar effects causing the modal changes of the distribu- 
tion function which can be interpreted as a nonequilibrium 
noise-induced phase transition. 

More complicated example is the rate function r(X) = 
a + bX and the input rate ([9]) . In this case 



w at (X) = Pq 



S(X) 



X 



-X/b 



exp{-fiX}(a + bX) 



A/6-1 



Pn 



1 + (fia/by x/b {X/b) exp{fia/b}r(X/b; fia/b) 
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Fig. 1. Stationary distribution function (I18|l in the storage 
model with r(X) — bX. Phase transition with increasing input 
intensity A. 



(r(x;y) is incomplete gamma- function). It combines two 
previously considered models: there is an atom at X = 
with weight Po and there is the phase transition at criti- 
cal A where the distribution switches from exponential to 
Gaussian-like. This critical value is A cr = a/i + 6. If 6 — > 0, 
it coincides with that for the model r = a, and if a — > - 
with the results of r = bX model 1|18|) . 

For a realistic release function r(X) = bX — cX 2 (l — 
X), (c, b > 0) corresponding to, for example, the nonlin- 
ear voltage-current characteristics, the solution to (|14| for 
exponential input ([9]) (expression l|15[)) yields [in this case, 
like for the linear model r = bX , Pq = 0]: 



g(x) = 



Nexp(-iiX)(X + c) x / b 



x exp 



{bX - cX 2 + cX 3 ) \cX 2 -cX + b | A /2fc 
Ac c(2X - 1) 



46 



arctan ■ 



a/c | c - 46 | 



, c < 46 



(we consider the case c < 46 only because natural restric- 
tions on the drift coefficient impose r(X) > 3J). The 
phase transition points can be explored from the custom 
analysis of the value of the corresponding cubic determi- 
nant Q. The number of phase transitions varies from one 
(Q > 0) to three at Q < depending on the relations be- 
tween 6/c, and A/6. This latter example can be com- 
pared to the nonequilibrium regimes found in the quar- 
tic potential well driven by an additive Levy-type noise 
[T4"I[1~5] ; as in the above systems the additional criticality 
is achieved due to a long-range character of the additive 
noise. 

The distribution functions and phase transitions in this 
class of models are common to various physical systems. 
Relative simplicity in mathematical treatment allows us to 
propose them as a handy tool for modelling physical phe- 
nomena of stochastic nature. We will show that this class 
of models based on the Poisson noise presents a prototype 
for stochastic modelling complementary to the commonly 
considered Gaussian noise. This latter, as the simplest case 



of phenomenologically introduced stochasticity, became in 
fact the most recognized way of introducing noise, and var- 
ious enhancements of stochastic description meant merely 
an extension of the Langevin source of a Gaussian nature 
introducing the small parameter of jumps value (see later 
in the Section 4). The use of the Poisson noise, starting 
from the storage model as its basic case, allows an exten- 
sion to more complicated cases as some regular develop- 
ment scries in a small parameter as well. In contrast to the 
Gaussian scheme, it does not require at the very beginning 
the smallness of jumps thus it is able of describing more 
adequately a wide class of physical phenomena where this 
assumption does not have its physical justification. As an 
example we mentioned a thermodynamic system of a small 
fizc, where the random value is the number of particles. 
The Gaussian assumptions are valid only for large size 
(in comparison to the rates of input or output) of a sys- 
tem, when the diffusion approximation can be used. The 
smaller is the system, the worse is the description in terms 
of minor random jumps (basic for Gaussian scheme), and, 
vice versa, more reliable becomes the description based on 
the Poisson character of a random process. 



3 Fluctuation-dissipation relations and 
formalism of the kinetic potential 

This section is a brief reminder of the formalism of the ki- 
netic potential [22| which is appropriate for presenting the 
properties of a Markovian random process in a compact 
form. 

The primary concept of the Markov process is the tran- 
sition probability for a random value B to jump from the 
initial state B\ at the time moment t\ into state B2 at 
that is the probability LO21 = w(i?2, t 2 \Bi, ti). The idea of 
Markovian behaviour imposes obvious restriction on the 
values W21, so that they possess a superposition property 
/ w 32 (B\B')u 21 {B'\B")dB' = w 3l (B\B"), (h < t 2 < t 3 ); 
in other words, form a continuous semigroup in time (the 
reverse element of this group for dissipative processes is 
not defined), so that all characteristics of a system can be 
derived from the infinitesimal generators of the group. In 
normal language, we bring into consideration the proba- 
bilities per time unit (l/r)u(B + A,t + r\B,t) (r -> 0). 
To characterize this function it is useful to consider its 
moments which are called "kinetic coefficients" : 



KJB, t) = lim - [ oj{B + A,t + t\B, t)A n dA . (19) 

r^0 r J 

The stationarity of the Markov process assumes that 
K n are time-independent. The kinetic equation for the 
distribution function of the process reads as [410.1,22 



du(B,t) 



()t 



E 

n=l 



d n \ K n {B) 



dB n J n! 



u(B,t); (20) 



w{B,t) = J uj ttl (B\B')Lj tl (B')dB' . 
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The kinetic potential is denned as the generating func- 
tion of the kinetic coefficients l22l: 



V(-6,B) = J2K n (B) 



{-or 



Thus the kinetic coefficients can be expressed as 
<:)' 



K n (B) = 



V(-9,B) {e=0 



(21) 



(22) 



for n = 1, 2, . . .. 
ten compactly: 



d{-6) 

With the equation $ZU§ can be writ- 



In ([23]) the notation Afd,B means the order of the differen- 
tiation operations: they should follow all actions with the 
multiplication by K n as it is seen from (|20l) . 

An example of the kinetic potential for the simplest 
and most utilized stochastic process is 



V{-9,B) = K l {B){-t 



\k 2 (B)9 2 



(24) 



With the choice K^B) = -b ■ B, K 2 {B) = D = const, 
K n> 2 = the corresponding kinetic equation is then the 
Fokker-Planck equation for the Orstein-Uhlenbeck process 
[4,11] which describes a system with linear relaxation to- 
wards the stationary solution in the Gaussian form 



y? 2 

~D 



uj s t(B) ~ cxp -— — 



Note that the kinetic potential for a process driven by a 
Levy flight noise has thus the generic form V(—9,B) = 
—QK\(B) + D9 a and assumes the kinetic equation of the 
formally fractional order which is not reduced to the series 
in ([20]): one uses instead its plausible generalization which 
can be encountered elsewhere (e.g., [T51[T71[T5 II 1 9 II2D] etc). 

As another example we write down the form of the 
kinetic potential for the class of storage models of Section 
2. The kinetic potential V(—9,B) through the transition 
probabilities of a Markov process is written as 



V(-9, B, t) = lim - \E ( e -s(x(t+r)-x(t)) | x , t ,\ _ 1 
t—>o t L V / 



E ( e -°l x *- x >) I X, 



^ x ^io tat2 (X 3 \X 2 )dX 3 
E(e- eX3 | X 2 )E(e ex *) 1 (25) 



where Xk = X{tk). Inserting there the Laplace transform 
of the random value from |T]), we obtain [7J, for an ele- 
mentary derivation see Appendix: 

V(-0,B) = -<p(0)+0r x (B), (26) 

where <p is defined in {7}, and r x {B) in (|3l4p . 

It is handy to introduce another generating function 
called the "image" of the kinetic potential [22]. Namely, 



let Wst(B) be the stationary solution u> = of the kinetic 
equation ([2"0"]) . The image of the kinetic potential V is 
defined as 



R(y,x) 



J exp{xB)w st (B)V(y,B)dB 



J exp(xB)uj st (B)dB 
or, in the notation of the transition probabilities, 



(27) 



R{y,x) = lim - x / / dB 1 dB 2 exp(xB 1 )(28) 

[exp(y(B 2 - B 1 )) - 1}uj{B 2 , t + T\B u t)w st {B x ) 
X / exp{xB)u st {B)dB 

The series of R(y, x) over y: 



Ey 
K n {x) — 
77.1 



n=l 



(29) 



defines new coefficients n n (x) being the image of K n : 



J K„(B)Lj st {B)exp(xB)dB 
J uj st (B) exp(xB)dB 



(30) 



K n (x) 



We note by passing that the variable x in (|27[) or (1301) . 
presenting merely a variable over which the Laplace trans- 
form of the process variables is performed, can be also un- 
derstood as a (fictive or real) thermodynamic force. This 
interpretation is clarified when we look at the "pseudo- 
distribution" exp(xB)uj s t (B) where xB stands for an amend- 
ment to the free energy of a system |22) . 

The reconstruction of a stochastic random process as- 
sumes that knowing macroscopic information about sys- 
tem we make plausible assumptions as to the fluctuat- 
ing terms of the kinetic equation, that is we try to con- 
struct the matrix of the transition probabilities ajy in any 
of equivalent representations (fl~9"]) , (|2"T]) or (j2"T[) leaning 
upon some macroscopic information about the random 
process. As the latter, we can understand the following 
two objects: 1) the stationary distribution ui st (B), and 
2) "macroscopic" equations of motion which are usually 
identified with the time evolution of the first momenta 
of uj(B,t) and hence the kinetic coefficient K\(B) (in the 
case of the sharp probability distribution where one can 
identify the "macroscopic variable" at all). For example, 
for the storage model scheme the problem is inverse to 
that considered in Section 2: knowing the macroscopic re- 
laxation law and the shape of the stationary distribution 
we then try to reconstruct the input function f(9). The re- 
laxation law is given by the balance of the averaged input 
and release rate p — r(X) (in the present class of storage 
models the input rate is X-independent, the generaliza- 
tion is considered further). Then, given r(X) one can set 
into correspondence to it the input function <p{9) yielding 
a given distribution uj st (X): 



V {9) 



J r(X)u s t(X)exp(-9X)dX 
J uj st (X)exp(-9X)dX 



(31) 
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The relations between said objects and the remaining 
part of the stochastic information contained in the process 
are called fluctuation-dissipation relations (FDR). These 
relations express the property of time reversibility of the 
transition probabilities (detailed balance). In the repre- 
sentation of the image of kinetic potential (|2T[> . (|28[) they 
are written in the most elegant fashion [22] : 

R(y + x, x) = R{— ey, ex) , (32) 

where e = ±1 according to the parity of the variable. The 
particular case of the FDR in the form ([32]) at y = 
represents the "stationary" FDR 

R(x,x)=0. (33) 

The FDR in the form ([33)) hold for any system in the 
stationary state with no assumption about the detailed 
balance, that is the system needs not to be in the equilib- 
rium state. Indeed, it is easy to check that ([33]) is just an 
another notation of the equation for the stationary distri- 
bution Mq,bV (-d/dB, B) uj s t(B) = (see Appendix). 

4 Reconstruction of the random process 

The problem of reconstructing a random process in the no- 
tations of the preceding section is formulated as a set of al- 
gebraic equations for a function R(y, x). Thus, given func- 
tions k±(x) and uj s t{B) we must find R(y, x) which identi- 
cally satisfies the relation ([32]) (for the system in equilib- 
rium) or (|33p (for the stationary system with no assump- 
tion about the thermal equilibrium and detailed balance). 
The problem, of course, has many solutions since those 
conditions do not define the function R uniquely. There 
exists the "FDR-indeterminable information" which hence 
should be borrowed from some additional criteria to be 
imposed on the equations in order to close the problem, 
which means confining ourselves within some class of the 
stochastic processes. The kinetic potential representation 
allows us to elucidate clearly the nature of the approxi- 
mations made. 

4.1 "Gaussian" scheme 

The standard reconstruction procedure considers the pos- 
sibility of setting the "FDR- indeterminable" functions neg- 
ligibly small, that is introducing a small parameter over 
which the kinetic potential can be developed in series [HI 
|2"2] . Thus the generalization of the "bare Gaussian" model 
showed in the example above ([24)) is achieved. In the series 

2 3 

V{y,B) = K 1 y + K 2 V -+K z V - + ... 

successive coefficients K n decrease progressively as 
This relation can be expressed introducing the family of 
kinetic potentials labelled by the large parameter /? [22] 
(compare with ([TT]) ') 

Vp(-0, B) = fiV(-6/P, B) , (34) 



This is a common approximation for a random process 
consisting in the fact that its jumps are small. If we keep 
only two terms K\ and K 2 , the second coefficient for one 
variable can be restored from FDR exactly. As an example, 
we write the kinetic potential reconstructed up to 4-th 
order (from the formula (|32[) applying development in the 
powers of y) : 

R g {y, x) = y Kl (x) (l - V - ) + y 2 (y - xf^- , 

where the coefficient k,±(x) is arbitrary (indeterminable 
from FDR). The first term describes the base variant cor- 
responding to (|24p. 

4.2 "Storage" scheme 

The assumption of small jumps leads to the possibility 
to neglect the higher order kinetic coefficients K n , con- 
structing a stochastic process by the "Gaussian" scheme 
(G-scheme). We suggest an alternative approach which 
can be regarded as complementary to the G-scheme and 
does not require the assumption of the small jumps. Like 
G-scheme, it has its basic variant which is well treatable 
mathematically. 

We assume now that the kinetic coefficients K n (B) are 
expandable into series over the variable B: 

B 2 B l 
K n (B) = k nfi + k Thl B + k rh 2 — + ■■■ k n,i-jY ••• • ( 35 ) 

Possibility of truncating these series implies that k n> i 
in (|35|) contain a small parameter 7 which decreases them 
progressively with the growth of the number I: k n ^i ~ 7' 
for n = 2, 3, 4, ....In the coefficient K\(B) determining the 
macroscopic evolution we however keep the macroscopic 
part r x (B) whose development on B does not depend on 
T- K^B) = -r x (B) +J2k 1>l B l /ll _ 

The image of the kinetic potential thus turns out to 
be a development into series 

R{y,x) = -y^(x)-Y jT[ (A l {x))<p l {y), (36) 
z=o 

here 

- T ^_ J r x (B)uj st (B)exp(xB)dB 
rx(X) J u st (B) exp(xB)dB 

, = fB l ij at (B)exp(xB)dB 
[ [X>> ~ Ju; st (B)exp(xB)dB ' 

and (pi are defined through coefficients in (|35[) : 
00 . 

-^(^E^- 1 = 0,1,2,... 

n=l 

The series ([3"6l can be considered as a development 
on the base {1, (A(x)), (A 2 (x)), . . .} which is a natural 
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base of the problem following from the peculiarities of its 
stationary distribution. The coefficient tpo(y) at y = — 9 
has the same meaning as the function Lp from ([7]) and 
truncating ([36| up to it reproduces the storage scheme 
of Sec. 2. This is a generalization referring to the multi- 
plicative noise processes dX(t) = —r x (X)dt + dA(t; B) in- 
stead of ([T]) with the characteristics function of the noise 
E(exp(-9A(t))) = exp(-t[ip + B(p 1 + .. .]) instead of ©. 

From the equation R(x, x) = applying it to (|36[) 
we obtain the reconstructed scheme for a random process 
("S-scheme"): 



Rs(y,x) 



OO ^ 

1=0 ' x 



(37) 



oo 1 

-y (r^xj- ^)) - £ -My)((A l (x)) - (A l (y))) . 
i=i 

Keeping in mind that the coefficient Ki(x) is supposed 
to be known the last expression can be rewritten as 



Ki{y) 
m{x) 



oo ^ 

^ T Mv){{A\x))-{A\ y ))) (38) 



The coefficients rji(y) at I — f,2,... are dissipative- 
indeterminable. They are given by 



(y) = -J2y n k n , l /n\ = vl (y)- yv > , l (o). 



which is now X-dependent (Sect. 4. 2). Find the solution of 
the simplest linear dependence of the input function on X 
with the kinetic potential (f2Tj) set as 



If we set ipi = for I = 1, 2, ... in (|36H38p , we recover 
the kinetic potential of the ordinary storage model (|26[) 
which can be restored from the stationary distribution 
exactly. The above formula (|31[) is a particular case of 
the application of this scheme. 



5 Conclusion. An example of possible 
application 

Both schemes sketched above - that is the common G- 
scheme and suggested complementary S-scheme of the stoch- 
astic reconstruction both lean upon two basic stochastic 
models which use respectively the assumption of Gaussian 
and Poissonian nature of the random noise. They both 
apply the series development of the kinetic potential on a 
small parameter. Keeping infinite scries leads to identical 
results in both cases. However, in real physical problems 
we use to truncate the expansion series keeping small finite 
number of terms. According to the physical situation and 
to the nature of the random noise either one, or another 
scheme would give a more reliable convergence. 

Consider now an example of the storage model with 
the linear release rate r = bX ([5]) and generalized input 



V(-6 ) X) = -(p (6) + 6bX-cXip 1 (6). (39) 



The first two terms in (f39|) describe the usual storage 
model with linear release, and the last term is the amend- 
ment to the input function proportional to X (cf. (|35l 
36)) ). The parameter c controls the intensity of this ad- 
ditional input. The equation for the Laplace transform 
F(9) = E(exp(—9X)) of the stationary distribution is 
V(-6,X ->■ d/d9)F(9) = 0, where the differentiation 
refers only to the function F. Its solution for V in (j3"9"| 



logF(0) 



ipo(u)du 
bu — c(fi(u) 



(40) 



For illustration specify now the input functions as 
^(6) = ^) = -^-, 



(41) 



which correspond to the exponential distribution functions 
of input jumps (|7l9p 6o,i( a; ) = /xexp(— fix). Then 



F(9) 



(i — cX/b 



-\/b 



(42) 



Comparing to the solution of the storage model c = 
shows that the additional term leads to an effective de- 
crease of the parameter /i — > /i — cX/b. The stationary 
probability distribution is given by the gamma-distribution 
function 



u st {X) = {^-cX/b) x / b X x / b - 1 e 



-l-X(n-c\/b) 



mx/p). 



The stationary solution exists if /j, — cX/b > 0, otherwise 
the system undergoes the phase transition with the system 
overflow likewise for the model with constant release rate. 
If c < 1 there are two phase transitions (increasing A), 
first of which is that of the model (fig)) (Fig.l) and the 
second is the system overflow at 1 — cp/b = 0. If c > 1 the 
overflow occurs earlier than the condition X = b meets and 
qualitatively the behaviour of the system is similar to that 
of constant release rate ((2]), with the transition condition 
cp/b = 1 instead of former p = a. 

Now let us sketch an example of application of the gen- 
eralized storage scheme for the problem of neutron fission 
process. Set the generating function ipi(9) in (f3"9"| in the 
following form 



Pi(0) = Ai 



fc=0 



7Tfc exp(— 6k) 



oo 
k=0 



here X)fc!Lo n kZ k is the generating function of the neu- 
tron number distribution per one elementary fission act 



Serge Shpyrko, V.V.Ryazanov: Stochastic Storage Models and Phase Transitions 



9 



(z = exp(—9) for a discrete variable); nk are probabilities 
of emerging k secondary neutrons (discrete analogue of 
the function b(x) in 0); Ai is the fission intensity (prob- 
ability of a fission act per time unit), Ai = (v)£f with 
average neutron velocity (v) and macroscopic fission cross 
section Sf [23]; further, b = l/l e f, lef is the average neu- 
tron lifetime till the absorbtion or escape [23]. Let us set 
c = 1 and assume that the function tpo (9) accounts for the 
external neutron source with intensity q = dtpo(9)/d9\g =0 
(the smallness parameter 7 in ()35|1 now describes the rela- 
tion A1/A0 of the intensities of fission and external source 
events). This probabilistic model is essentially the same 
as in the example of the generalized storage scheme (|39|) 
sketched above. From (l23l) , (l39l) 



dlogF(0)/dt = {-MO) - [9b-^(9)}d/d9}\ogF(9) , 

that is we arrive at the equation for the distribution func- 
tion of the prompt neutrons in the diffusion single- velocity 
approximation. The macroscopic equation for the averages 
(N) = -dlog F(9)/d9 w=0 is 



d(N) 
dt 



[(v){v)E f -l/l ef ] (N)+q 



|23| and coincides with that obtainable apriori from the 
stochastic storage model. The neutron reproduction fac- 
tor is defined as k = Xi(v)/b = {v){v)Sfl e t [53], where 
(i>) = \^ 1 difi(9)/d9\g =0 is the average number of sec- 
ondary neutrons per one fission act. The expression for the 
generalized storage model phase transition \i(v)/b = 1 
corresponds to the reactor criticality condition k = 1. 
Extending probabilistic schemes in (l39|) beyond the toy 
model considered here and introducing vector (multi-comp- 
onent) stochastic processes allows for taking into account 
the delayed neutrons, as well as various feedbacks and 
controlling mechanisms. 



6 Appendix. Derivation of the relations (125] ) 
and (S3]) 

Here we sketch the elementary derivations of the expres- 
sions in the text describing the kinetic potential of the 
storage model (expr. ([2"rj]0 and that of the stationary FDR 
through the image of the kinetic potential (expr. ([33]) ). The 
more rigorous and generalizing derivations can be searched 
elsewhere, resp. 0|3l[7l[22] . 

6.1 Kinetic potential of the storage model 

From (0 and ([25]) using the fact that the input rate is a 
random process independent on X(t) (to simplify nota- 
tions we set the initial moment t — 0): 



f 



-E e 



-SA(r)+e ^r x [X(u)]du _ j 



i S.E (e- 0A ^ E (l + 9 £ r x [X(u)]duj - 1 j 



Then, using (O and taking t — > 0, 

V{-6,B) = -<p{6)+ 9 lim -( / r x [X(u)]du) . 

The last term of this expression gives #r x [AT] + 0(t) in 
the limit r — > provided that the intensity of jumps is 
finite which is the case of the considered class of Poisson 
processes. 

6.2 Stationary FDR 

We limit ourselves to the noncquilibrium FDR relation 
only. For a general case of the detailed balance, as well as 
for the generalization to non-Markov processes, the reader 
is referred, e.g., to the book [22] , 

The stationary Fokker-Planck equation is written as 



51 n\ 



dB 



K n {B)u st {B) = 
(43) 

Perform over f4"3"]) the operation J exp(xB)(-)dB. Use the 
relation 



_d_ 

~dB 



f(B)dB 



e xB f(B)dB 



for some f(B) which can be verified with recursive inte- 
grations by parts (the terms with full derivation vanish 
at B — ±oo; if the space of states is a semiaxis as in the 
storage models, the integration J Q _ including the atom in 
is assumed). Then, 



O j e xB cu st (B)dB (f^ ~ } x n K n (B) 



R(x,x) = 0. 
(44) 
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